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The evolution of a scalar field is explored taking into account the presence of a background fluid in 
a positively curved Universe in the framework of loop quantum cosmology. Though the mechanism 
that provides the initial conditions for inflation extensively studied in the literature, is still available 
in this setup, it demands that the initial kinetic energy of the field be comparable to the energy 
density of the background fluid if the fleld is initially situated at the minimum of the potential. It is 
found, however, that for potentials with a minimum such as the chaotic inflation model, there is an 
additional mechanism that can provide the correct initial conditions for successful inflation even if 
initially the kinetic energy of the field is subdominant by many orders of magnitude. In this latter 
mechanism the field switches direction when the Universe is still in the expanding phase. The kinetic 
energy gained while the field rolls down the potential is subsequently enhanced when the universe 
enters the collapsing phase pushing the field one step up the potential. This behavior is repeated 
on every cycle of contraction and expansion of the Universe until the field becomes dominant and 
infiation follows. 



I. INTRODUCTION 



Currently, the leading background independent and 
non-perturbative candidate for a quantum theory of grav- 
ity is loop quantum gravity P, _3J which is a canoni- 
cal quantization of general relativity based in Ashtekar's 
variables. This approach provides a discrete structure 
of geometry. Continuous spacetime emerges in the large 
eigenvalue limit of quantum geometry. Loop quantum 
Cosmology (LQC) is the application of loop quantum 
gravity to homogeneous and isotropic mini-superspaces 
An important feature of LQC is that eigenvalues of 
the inverse scale factor operator are proportional to pos- 
itive powers of the scale factor below a critical scale a* 
In particular, it was shown that this property af- 
fects the behavior of the kinetic energy of a scalar field 
at small scales. There has been considerable interest in 
understanding the dynamics of a scalar field in this semi- 
classicalphase, in particular, its role in the origin of in- 
flation IE II S, H ES El 111 lilJ3' avoidance of a 
big crunch in a closed Univers e Il5l Il6| and non-singular 
bounces in the cyclic scenario [l7|. 

In classical standard cosmology the energy density of 
a perfect fluid with constant equation of state evolves as 
p oc a-3(ui<,i+i) i^Qnce , diverging as a ^ for Wd > —1. 
Consequently, the presence of a matter component into 
the dynamics may lead to singularities. The behavior of 
matter in the semi-classical phase, however, is not com- 
pletely understood as a theory of quantum gravity includ- 
ing matter is yet to be constructed. A purely phenomeno- 
logical approach has recently been attempted Ell- It 
was found that similarly to the case of a scalar field, the 
classical cosmology equivalent of the energy density of a 
perfect fluid is modified in LQC. It varies proportionally 



to positive powers of the scale factor (for Wc\ > 1/4) at 
small scales, preventing it from diverging. 

The semi-classical phase of LQC is defined for values 
of the scale factor in the interval < a < a*, where 



and al = jaf/3. The Barbero-Immirzi pa- 
rameter is 7 = ln2/\/37r « 0.13 and j is a half in- 
teger quantization parameter. The evolution equations 
are modified by the presence of non-perturbative effects. 
The discrete nature of spacetime is important below a.; 
and is described by a difference equation. The evolution 
equations take the standard classical form above a*. 

The modified Friedmann equation for a positively 
curved Friedmann-Robertson- Walker Universe sourced 
by a scalar field (p with self-interaction potential V{(j)) 
and a background fluid with classical equation of state 
Wci, has the form 



= y (P0 + Ph) 



1 



(1) 



where we have defined k = Snlp^ — SttTOpj and the 
energy densities of the scalar field and background fluid 
as El 
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(2) 

Pb = PoU^"a _ (3) 

The quantum correction function D{q) is defined by @ 



D{q) = (-j g^/^{7[(, + l)"/^-k-ll"/^' 



llq 



(g+l)V4_sgn(g-l)|9-l|V4]|'(4) 
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with q = (a/a*)2. In the semi-classical phase (a <^ a*), 
the quantum correction function varies as D oc a^^ , has 
a maximum near a — and falls monotonically to the 
classical limit 13 = 1 in the classical phase (a > a*). 



2 



The equation of motion for the scalar field is given by 
1 dlnD' 



3H 1 



3 din a 



(5) 



Diff'erentiating Eq. and substituting for cj) using 
Eq. lO gives 



H = — 



1 - 



1 dlnD 



K 

Y 



6 din a 

Wci d\nD 



3 din a 



(6) 



The effective equations of state of the scalar field and 
background fluid can be written as 0| 
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02 + 2DV 
^ 1 dluD 
3 d In a 



1 dhiD 
6 din a 



(7) 
(8) 



We can see from these equations that deep into the semi 
classical phase, the equations of state become Wff, ~ —4 
(neglecting the potential) and Wh = —'iwd- It was shown 
in Ref. il5j| that such a feature allows the energy density 
of the scalar field to cancel out the curvature term al- 
lowing the Universe to re-bounce at small values of the 
scale factor, hence avoiding a big crunch. A recoUapse 
also occurs at a large value of the scale factor in the clas- 
sical phase as the kinetic energy of the field scales more 
quickly than the curvature term. An oscillatory Uni- 
verse is therefore the natural outcome of the LQC cor- 
rections in a positively curved space sourced by a scalar 
field. Since the kinetic energy of the field never van- 
ishes while the potential is negligible, it was discussed 
in Refs. UlIS Q 113, that the oscillatory behavior of 
the Universe enables the scalar field to move up in its 
potential establishing the initial conditions for slow roll 
inflation (when cjp k, V), even admitting that the field 
starts its evolution at the minimum of the potential with 
small kinetic energy. 

In this paper we extend these works by including a 
background fluid in the dynamics of the Universe. We 
find that a mechanism similar to the one just described 
is still a possibility but it only results into an inflation- 
ary expansion for a sufficiently large contribution of the 
field's kinetic energy. Nonetheless, we describe an alter- 
native mechanism through which the field can become 
dominant, even if the kinetic energy is initially many or- 
ders of magnitude below the energy density of the back- 
ground fluid, and discuss the parameter space where it is 
most efficient. 

Under the notion of perfect fluid, it is typically as- 
sumed that quantum effects play a negligible role in the 
dynamics of the corresponding quantum field such that 
the concept of an equation of state can be applied. Since 
in our study we are precisely concerned with evaluating 
the effects of loop quantum corrections on a fluid with a 



constant classical equation of state, we are in danger of 
dealing with an inappropriate definition. We must nec- 
essarily assume that the dynamics affected by the loop 
quantum modifications has a time scale much larger than 
the time neces sary to ensure thermodynamical equilib- 
rium (see also [iMEl)- This assumption, however, may 
become invalid when the scale factor of the Universe be- 
comes close to ai where the full quantum theory comes 
into play. We will see, however, that the interesting and 
viable models including a background fluid are only well 
defined for a» ^ ai and a « a*, hence, we expect that 
the quantum corrections exert a harmless role where the 
definition of equation of state and energy density is con- 
cerned. 



II. CRITICAL POINTS AND STABILITY 

The equations of motion |SJ| and © can be rewritten 
as a system of first order differential equations by defining 
the dimensionless quantities 



V3 



\V\a, 



z 
w 



(9) 
(10) 

(11) 
(12) 



For generality, we have considered the possibility that the 
scalar potential and/or the energy density of the back- 
ground fluid are negative (e.g. a negative cosmological 
constant). The system is now governed by the following 
equations: 




1 ± 



D \xy , 



- 1 ± 



±^u;2 [u;^j(A_3)_i] 



2qz, 

-VODX^xiT-l), 
subject to the Friedmann constraint 



9 
A' 



2 I 2 

X ±y 



2 I 2 
Z ± W 



1 



(13) 
(14) 

(15) 
(16) 
(17) 

(18) 



Here a prime corresponds to differentiation with respect 
to conformal time dr = dt/a. We have also used the 
definitions: 



A 



1 j_ dy 

kVcM 



(19) 



3 



d^V fdV 
dlnD 



r = V I — 
A = 2 



din g 

It is useful to define the dimensionless quantities 



r = y Po a. 



-1 — SlWcl 



V = —Vai 



(20) 
(21) 

(22) 



so that y and w can now be written as 



±w^ = rD"""' q-(i+3««ci)/2 _ (23) 



These definitions will allow us to perform a general study 
of the system regardless of the value of the quantization 
parameter j which is now enclosed in r and v. 

We will assume for now that the potential is constant, 
i.e. A = 0. In the next section we will extrapolate the 
results to a varying potential under the assumption that 
the variation is slow enough such that we can refer to 
"instantaneous critical points". Following IJj we are 
looking for static solutions such that d = a = or, in 
the current variables, q' — and z' = 0. The system has, 
therefore, critical points in Zc = and qc such that 

± 7/2(6 - A,) ± wl (Ae - 3){wci - 1) + A, - 4 = . (24) 

Note that Wc is the value of w defined in Eq. (|12|l at the 
critical point and Wd is the classical value of the equation 
of state. 

We must also guaranty that the kinetic energy of 
the scalar field is non-negative at these critical points, 
hence it follows, using the Friedmann constrain H18|l and 
Eq. ^ that 



1 



6- A, 



[2±wl{w,i{A,-3)-3)]>0. (25) 



In Fig.Qwe show, by dotted lines, the relation between 
V and q of the system's static solutions as determined 
form Eq. (|24|l . The case studied in Ref. corresponds 
to the top lines, r = 0. 

The nature of the critical points can now be determined 
by linearizing q' and z' around the critical points. The 
eigenvalues of the system are given by 

= gc^(9c) [lTyc±Ki-l)u;c] ±2/c(6- A,) 

T^wl Ki(Ac - 3) - 1] (A, - 3)(wei - 1) . (26) 

Stable solutions exist when < 0. Substituting for 
yc in Eq. H26|l using Eq. H24|l we observe that the eigenval- 
ues are regular everywhere except at qc = 0.835 (where 
A = 6) and qc ~ I. These are the points where yc and 
dA/dg blow up and change sign, respectively. Therefore, 
a transition from a stable critical point (m^ < 0) to a 
saddle point (m^ > 0) is expected at these values of qc- 
Of course, transitions may occur at additional points. 



2.5r 
2 

1.5 
1 

0.5 


-0.5 

-1 - 



0.5 




1.5 



FIG. 1: Dependence of v on the critical value qc- From top to 
bottom the curves represent r = 0, 0.3, 0.6, 0.9 and 1.2. Solid 
lines represent stable critical points and dashed lines represent 
saddle points. Dotted lines correspond to the position of the 
critical points if condition <\25^ was satisfied. 



We depict these transitions in Fig. ^ by representing the 
stable points with solid lines and the saddle points with 
dashed lines. The classical equation of state is assumed 
Wci = 1/3. We will always use this equation of state 
when dealing with specific examples in the figures, but 
keep the analysis general in the analytical derivations. 



III. PHASE SPACE TRAJECTORIES 

In this section we will turn our attention to evaluate 
the trajectories of the phase space. The first point to keep 
in mind is that the kinetic energy of the field must always 
be definite positive, i.e. x'^ > 0. Using the Friedmann 
constrain we can rewrite this condition as > zLy^ ± 
— 1. On the other hand, since 2;^ > we have that 
in the regions of the phase space where zty^ ± — 1 < 
the latter condition on is automatically satisfied. 
Conversely, in the regions where zty^ ± — 1 > 0, there 
are regions of exclusion for Ha^ ~ z/^, with upper and 
lower bounds defined as 



b < Ha^ < b , 



(27) 



where & = {±y'^±w'^ — iy^^/y/q. By inspecting these con- 
straints we can have an idea of how many regions of ex- 
clusion to expect given a pair (v,r). First we recall from 
Eq. (j2Sll that increases linearly with q. We note that 
Ph oc g6M)ci-3/2 diverges as q ^ when Wd < 1/4, hence, 
we will not consider these cases from now on. Since, 
the semi classical phase it follows that 
increases with q in the semi classical phase. More- 
over, w'^ decays away as w"^ oc g-(i+3™ci)/2 for g > 1. In 
general there are four cases to consider. They correspond 



FIG. 2: Areas of the parameter space {v, r) where there exists 
one region of exclusion (unshaded area) and where there are 
two exclusion regions (shaded area). Regions of exclusion are 
regions in the phase space where the kinetic energy of the field 
is negative. 



FIG. 3: Illustrating the phase space trajectories (solid lines) 
and regions of exclusion (shaded areas) for r = 0.6 and vary- 
ing V. The trajectories evolve counter clockwise. 



to the various combinations of signs before and w'^ in 
the definition of b. We represent those combinations by 
the pair (i, j) where the first element represents the sign 
before and the second element, the sign before w^. 

(i) When we have the pair (+, +), we can have up to 
two regions of exclusion which merge as either w or r 
increase. One region appears at small but non-vanishing 
q and the second extends from small q to infinity. 

(ii) When we have (+,—), again we can have up to 
two regions. This time, they merge as the ratio \v/r\ 
increases. The regions are qualitatively like the ones in 
(i). 

(iii) When the choice of signs is (—,+), there is only 
one possible region of exclusion which exists for small 
non-vanishing q. The area of the region increases with 
increasing ratio \r/v\. 

(iv) Finally for the combination (—,—), there is no 
region of exclusion and the whole phase space is available. 

This is a qualitatively description that can be comple- 
mented by a quantitative analysis shown in Fig. [21 The 
shaded area corresponds to choices of {v, r) that lead to 
two regions of exclusion in the phase space of the system. 

In FigO we illustrate how the areas of exclusion and 
trajectories evolve if we fix the value of r and increase 
the potential. Indeed, Fig. corroborates the informa- 
tion given in Fig. |21 More specifically, setting r = 0.6, 
we verify that up to w « 0.31 there is only one region of 
exclusion, however, a second region appears above this 
value. This region increases in area as v is further in- 
creased and merges with the first region for v w 0.4. 

Let us now look at a specific example in the context 
of a realistic varying potential. We assume that the field 



is initially (/)init = and does not go any further than 
(j) = 3mpi. Again we choose r — 0.6. From Fig. ^ we 
have that for such a value of r, the maximum value of 
V allowed by the kinetic constraint is Wmax = 0.31. For 
a quadratic potential V = m^0^/2 with m ~ lO^^mpi, 
we have in order to concretize the point of maximum dis- 
placement, (p = 3mpi, that the quantization parameter 
must be j = 1.94 x 10^^. Moreover, we can also read from 
Fig. Glthat when V = 0, qc ^ 0.866. Choosing this value 
for ginit along with the other parameters, we obtain the 
evolution shown in Fig. 01 Indeed we see that the poten- 
tial and the scale factor are oscillating in two fashions. 
One around the instantaneous critical point and the sec- 
ond between the critical points corresponding to V — 
and V = Svmax/ K^oi- We can understand this behavior 
by returning to Fig. |21 The evolution starts at u = 
and develops along the line r = 0.6. As the potential in- 
creases, a second region of exclusion is generated at the 
position of the critical point when v = 0.31. This re- 
gion increases as the value of the potential increases. In 
fact, it can become as large as the orbit of the trajectory, 
in which case, by the definition of the regions of exclu- 
sion, x'^ = 0. This means that the field has reached the 
point of maximum displacement and turns around. The 
potential decreases and so does the area of the bounded 
region of exclusion. Once the field passes through zero 
and changes sign, the evolution mimics the evolution on 
the opposite side of the potential i.e. the field enters a 
cyclic regime. 
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FIG. 4: Evolution of q and the potential V, with time. The 
initial conditions are i/iinit = 0, //init = 0, qi„it = 0.866, j — 
1.94 X 10" and r = 0.6. The V axis is labeled in Planck units. 



IV. CONSEQUENCES FOR INFLATION AND 
"GRACEFUL ENTRANCE" 

We have seen in the previous section that, in some 
cases, as the potential increases, an exclusion region ap- 
pears in the phase space which can cause the field to 
loose all its kinetic energy and therefore to turn around 
in the potential, however, in a regime of static solutions 
d = a = 0, not in an accelerated expansion, d > 0, as 
required for inflation. Looking at Fig. El we see that pro- 
vided < r < 0.38, no second region is formed. The 
upper bound corresponds to the value of r above which 
x^{q — 1) in Eq. (|25l) changes sign, or in other words, is 
the lowest value of r for which the corresponding solid line 
and dashed line (like the ones in Fig. ^) do not touch at 
g = 1. Hence, the area of parameter space < r < 0.38, 
allows in principle, a smooth transition between a regime 
where the field is being pushed up the potential to a phase 
where this solution becomes unstable and the field slow 
rolls back down the potential thus making the Universe 
to inflate in the very same spirit of Refs. |9llT^ll^. 

In this article we are primarily interested in those situ- 
ations in which the background fluid is initially dominat- 
ing the evolution of the Universe, and the curvature term 
is important only at the values of the scale factor where 
the bounces occur. More specifically, neglecting the con- 
tribution of the scalar field, we require at the bounces 

Ph{a = ai^2)/3 ~ aj~2 (or equivalently ±u'2.2 = 1)1 
which leads to the constraint 



' ^1,2 yi,2 



= 1 



(28) 



This equation has two real roots (qi and 92) provided 
r > 0.82. It is important to stress that this value lies 
outside the bound < r < 0.38 that avoids the genera- 
tion of a second exclusion region. This means that it is 
not trivial to implement a mechanism akin to the one dis- 



cussed in Refs. 0, \lAJ when a dominant background 
fluid is present. In fact, we can only expect a transition 
from an oscillating Universe into an infiationary one by 
admitting that the amplitudes of the oscillations of the 
scale factor are sufficiently large such that the trajecto- 
ries in the phase space touch the separatrix (where the 
trajectories appear to intersect) before the diameter of 
the second exclusion region equals the diameter of the 
trajectory. 

What is therefore the correct values of initial param- 
eters that deliver the desired evolution? Admitting that 
initially the evolution starts with V — and H — 0, we 
have xf^^^ + wf^^^ = 1. On the other hand, when the field 
turns around in the potential to slow roll back down in 
the last cycle, very near when the two exclusion regions 



merge, one has y. 



2 

merge 



W„ 



1. Further admitting 



that qi = ginit ~ 9finai and defining the useful quantity 



-Dinit Pbinit 



(29) 



that relates the kinetic energy of the field to the energy 
density of the background, using Eq. H28|l and requiring 



>yn 



we find that 



a>2 ginitf n 



(30) 



It is concluded that the ratio of kinetic energy to the 
background's energy density is bound from below by a 
quantity that is independent of j . More specifically, we 
realize that it increases for increasing qinit- Indeed, from 
Eq. when qi (= ginit) increases, r decreases, and 
Fig. [51 implies that Emerge increases. For example, if 
^init = 4 X 10"'', the two region of exclusion merge for 
u w 6 X 10"^ which results into the limit a > 10"^^°. It 
is however worth pointing out that if qi^n is small the 
consistency condition Hai < 1 is violated as can be seen 
from Fig.ini This constraint simply states that the Hub- 
ble radius must be larger than the limiting scale of the 
theory. On the other hand, for large ginit, the ratio a 
becomes close to unity meaning that the contribution of 
the field must be important. 

In the remainder of this paper we deal with an alter- 
native mechanism where the energy density of the scalar 
field, though initially subdominant, builds its way up 
during the numerous cycles of oscillation of the Universe 
until it becomes dominant delivering an accelerated ex- 
pansion. We show in Fig.[5lan example of such an evo- 
lution. The example yields only a few e-folds of infla- 
tion, for illustrative purposes, but as we shall see below, 
there are regions of parameter space that provide a suc- 
cessful inflationary scenario with more than 60 e-folds 
of expansion. The evolution of the field starts in the 
semi-classical phase, Inginit = —7.8, at the bottom of the 
potential with a kinetic energy that is twenty orders of 
magnitude smaller than the energy density of the back- 
ground fluid. The anti-frictional term in the equation 
of motion of the scalar field accelerates the field up the 
potential as the Universe expands. The energy density 
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FIG. 5: Evolution of the energy densities with the scale fac- 
tor. The scalar field is represented by the solid line, the back- 
ground fluid by the dashed line and the curvature term by 
the dash-dotted line. We have used a — 10"^", <^init = 0, 
Qinit = 4 X 10~* and j = 9 x 10^ for a quadratic potential 
V = rr?(^ jl with m = 10~^mpi. The vertical axis is labeled 
in Planck units. 



At this point we need to establish which set of the pa- 
rameter space leads to a viable model of inflation rather 
than to an infinite series of oscillations of the Universe. 
Let us start by calculating by how much the scalar field is 
displaced in half a cycle of the Universe. We assume that 
the initial value of the scale factor is deep into the semi- 
classical phase and that the potential is initially negli- 
gible. From the equation of motion of the scalar field, 
neglecting the contribution of the potential, we obtain 



( ^ 



-J-^ — V^init 

d In a V flinit 



D 



A, 



(31) 



Using iP' « K^pb/Z (this is obviously not true at the 
bounces but it is a good approximation within the ex- 
tremes of the scale factor) and integrating one gets 



(^/ = 0init + 



2 gjnit"* ^mit 

3 AnitV^ (9-4ii;cl)(l -Wcl) 



X <^ (10 - SWcl) 



3(l-t«ci)/5 



of the field evolves with an effective equation of state 
Wcf, = —4. Once in the classical phase (a > a*) this 
term turns into a frictional component forcing the field to 
slow down. The equation of state, therefore, approaches 
— 1. In the case of Fig. [S] the field's kinetic en- 
ergy decreases below the potential energy at In g/ « 7.8 
and comes to a stand and turns around in the potential 
{wtf, = —1) when In^f sa 9.7 (indexes "/" and "i" stand 
for "freeze" and "turn around", respectively). On the 
way down the potential, the field gains again kinetic en- 
ergy which is enhanced once the Universe recollapses (as 
in the classical phase the frictional term becomes anti- 
frictional when the Hubble ratio is negative). The field 
is again slowed down as soon as the scale factor enters 
the semi-classical phase during the collapse. When the 
kinetic energy becomes smaller than the potential, the 
equation of state again approaches —1. In this example, 
the kinetic energy decreases but it does not vanishes in 
the semi-classical phase; hence, when the Universe goes 
through the bounce starting a new cycle, the field keeps 
moving in the same direction. The subsequent evolu- 
tion follows closely the previous description for a number 
of oscillations of the Universe. The remarkable feature 
about the evolution is that the field gains energy on each 
cycle thus allowing it to be displaced by a larger amount 
from one cycle to the next. When the contribution of 
the energy density of the field becomes important, the 
field can slow roll and the Universe naturally starts an 
inflationary evolution say, in a " graceful entrance" . 



-(1-Wci) f y 



3(l-u)ci) 



(27-12u;,i)/4 
Vinit 



(32) 



where (j)f and qf correspond to the value of the field and 
q when the field freezes in the potential, respectively. To 
perform this integration we have considered two epochs of 
evolution. In the semi-classical phase, between ajnit and 
as the function D takes the asymptotic form -Dapprox = 
(12/7)^ g^^/^. In the classical phase, when a > 05 we 
take D — 1. The quantity ag is defined as the value of 



the scale factor where D„ 



1. We can estimate the 



value oi qf by using the expression 



0/ = 0init 



6 3/2 
<it 9/ 



(33) 



which follows from using the first integral of the equa- 
tion of motion of the field along with the asymptotic 
Furthermore, a freezing point can only 



form D 



approx- 



be identified provided the potential becomes at least as 
large as the kinetic energy, such that, m'^(j)j/2 « (pj, for 
a quadratic potential. Using Eq. (|33|l we obtain 



r;3 



7 



2 ( _ 1 'Punt 
12 



(34) 



Using Eq. (|29|) we can write the initial kinetic energy as 



3 a r 



-3(l+u.„i)/2 jjl+w,i 



(35) 
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FIG. 6: The dashed hne represents the lower hmit on j result- 
ing from requiring that the field freezes before the bounce at 
§2. This is estimated using Eqs. 13211 and 13311 . The solid line 
represents the borderline between an evolution that pushes 
the field up the potential in successive steps (below the curve) 
and an evolution where the field is displaced from the initial 
position stops before the bounce, turns around and is dis- 
placed again in the opposite direction in the collapsing phase 
(above the curve). This curve was obtained numerically and 
fits to g/ = g2'*'^ accurately. The dash-dotted line gives a 
lower bound on the values of j that satisfy the consistency 
condition Hai < 1. 



Inserting this expression into Eq. (|32|l . the multiplicative 
factor before the curly brackets reads, 



(to,i-l)/2 3(l-«)d)/4 
Vinit 



(36) 



which is explicitly independent of a*. Moreover, when 
/2 the qf dependent term in Eq. I|32|) 



> 



can be neglected and then the displacement (pf is com- 
pletely independent of the value of j for a given qinn and 
a. From Eq. (|35|l and substituting for r using Eq. H22|) 

we get oc a* hence we expect (^/ to become 

independent of j for small values of this parameter. 

A necessary condition for the field to turn around is 
that Qf < (?2, i-e., it must freeze before the Universe rec- 
oUapses. This constraint results into a lower bound on 
the value of the quantization parameter j with respect to 
the initial value of ginit- We draw this bound in Fig. El 
dashed line. This condition might not be sufficient to 
guaranty that the field stops and reverses its direction, 
as the kinetic energy can decrease below the value of the 
potential but without vanishing. In this latter case the 
field keeps moving up the potential in a series of steps 
before the turn around occurs (such cases were discussed 
in the first part of Section HVI an example is illustrated 
in Fig. ^ and do not result into an infiationary expan- 
sion, as we have seen). Understandably, Ingy must be a 
fraction of the value of In (72- By numerically integrating 



the equations of motion it can be verified that this frac- 
tion is nearly a constant and its value is 0.82. In Fig|H| 
we show by a solid line the lower limit on j for which the 
field does changes direction on each cycle corresponding 
to Qf < 

We also draw in Fig. |S1 the limit on j imposed by the 
consistency condition Hai < 1. Using the approximation 
that the maximum of the background energy density oc- 
curs when q = I and D = I, this condition is equivalent 
to j > 3 r. Using qi = qiait in Eq. H28|) and the asymp- 
totic form D 



approx 5 



we get that deep in the semi-classical 
phase the limit on j can be explicitly written in terms of 
the initial value ginit as 

(37) 



(l-12t«d)/2 
yinit 



From Fig. 1^1 it is concluded that below ginit = 2 x 10^'' 
the consistency constraint establishes the most important 
lower limit on j. 

Finally we can evaluate how efficient the mechanism 
is. By numerically integrating the equations of motion 
it is possible to determine how high can the field move 
in the potential before it rules the dynamics of the Uni- 
verse. More specifically, we have evolved the system up 
to qcnd = 1 + (In 52)72, i.e. until the scale factor is one e- 
fold larger than its maximum value at the bounces, when 
the background fluid is dominant. We have assumed a 
quadratic potential V = m'^0^/2 with m = 10~^mpi. 
Figures [T] and |S1 show the maximum value of the field 
and the maximum value of the Hubble ratio, respectively, 
with respect to the quantization parameter j. In all cases, 
the evolution starts from 0init = 0. 

Important information can immediately be extracted 
from these figures. First, we note from Fig. |S1 that the 
consistency condition Hai < 1 (which is equivalent to 
< 7.7 in Planck units) is only satisfied provided 
j > 8 X 10'', j > 3 X 10'' for the first and second ex- 
amples illustrated. Above these limits there are values of 
j for which the maximum value of the field is above three 
Planck units, hence providing a viable inflationary sce- 
nario with more than 60 e-folds of accelerated expansion. 
We note, however, that by increasing gi„it the maximum 
achieved value of (p decreases, as it is seen by comparing 
the middle and bottom panels of Fig. ITJ. 

A remarkable feature about these figures is that they 
present a tooth wood saw bladelikc structure for low j 
values. This can easily be understood qualitatively. In- 
deed, from Eqs. and substituting for r using 
Eq. (|22|) . we see that the value of qf decreases for increas- 
ing j as qf oc a;7^ (under the approximation that (j)f 
does not depend on j, which is accurate for (pf^-^^^ ^ (pj 
as we have seen before). This means that in a Universe 
where j is large, the field freezes and consequently starts 
rolling down when the scale factor is smaller than in a 
Universe with low j. The field is accelerated during the 
expansion phase between qt and 52 (even though the 4> 
term in the equation of motion is a frictional one) , dur- 
ing the collapse between 52 and q ~ 1, and again during 
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FIG. 7: Maximum achieved value of the field (in Planck 
units) with respect to the quantization parameter j for three 
different sets of initial conditions. Upper panel: ginit = 4 x 

ginit = 4 X 10~* and 
= 2.5 X 10"^ and a = 
i = 3mpi. The vertical 



10^ and a — 10~ ; middle panel 
a — 2 X I0~^^; lower panel: ginit 
2 X I0~^^. Dashed lines represent 
axes are labeled in Planck units. 
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FIG. 8: Maximum achieved value of the Hubble ratio (in 
Planck units) with respect to the quantization parameter j for 
three different sets of initial conditions. Upper panel: ginit ~ 
4 X 10"* and a = 10"^°; middle panel: ginit = 4 x 10~* 
and a = 2 X 10~^^; lower panel: ginit = 2.5 x 10~^ and a — 
2 X 10~^^. Dashed lines represent the consistency constraint 
Hui < 1. The vertical axes are labeled in Planck units. 



the expansionary phase when q < 1. But since the field 
starts moving ear her in a Universe for larger j, the net 
effect is that it is pushed further up the potential in this 
case. In other words, the field gains kinetic energy faster 
if j is large, which results into a larger displacement of 
the field from its initial position. This explains why the 
maximum value of the field increases with increasing j in 
Fig. The reason why the maximum value of the field 
suddenly drops by increasing j is understood by realizing 



that for a critical value jc there is one less cycle than if 
j ^ jc, thus preventing the field to move as high in the 
potential. As we keep increasing j, the description made 
above follows, now with one less cycle of expansion and 
recollapse. 

Along the same lines, a similar oscillatory structure is 
expected by varying a. We have seen in Fig [S] that for a 
given j the energy density of the field increases on each 
cycle by a nearly constant amount. It can be verified that 
if a was instead 10~^^ (10"^'*) the evolution would nearly 
reproduce the one in Fig. |21 with the only exception that 
there would be one less (more) cycle. 

Another prominent feature in the figures are the sharp 
peaks giving place to smooth extremes above some value 
of j. Indeed, when j becomes sufficiently large, the field 
has enough time to increase its speed before the recol- 
lapse takes place. We have in those cases (jp > V ai 
q = q2. Therefore, the field's equation of state is consid- 
erably larger than = — 1 at this point. This behav- 
ior translates into a similar effect of the one of having a 
slightly smaller j, i.e. the ratio of energy densities from 
one cycle to the next is smaller than expected had the po- 
tential energy remained always dominant. As j becomes 
even larger, the field is able to perform a few oscillations 
around the value at the minimum of the potential before 
the collapse occurs. To summarize, the effect of increas- 
ing j is compensated by the fact that the field loses most 
of its potential energy, thus the number of cycles of the 
Universe before the field becomes dominant does not de- 
crease anymore as j is increased. 



VI. DISCUSSION 

In this article we have looked upon the dynamics of a 
scalar field evolving in a closed Universe where a back- 
ground fluid is present in the context of loop quantum 
cosmology. In the first part we have searched for static 
solutions of the system and evaluated their stability ad- 
mitting a constant potential. We have discovered that, 
for a constant value of po and j, as the potential is in- 
creased, a second region of exclusion (where the kinetic 
energy is negative) is generated in the phase space. We 
argued that the occurrence of such a region makes it diffi- 
cult to implement a similar mechanism to the one studied 
in Refs. where the field moves up the potential 

in a series of steps resulting from the cyclic expanding 
and contracting phases of the Universe. A viable model 
requires the phase space trajectories of the system to be 
sufficiently away from the critical point in order to avoid 
the exclusion region since its area grows as the potential 
increases. This requirement translates into a lower limit 
on the initial value of the kinetic energy of the field and 
is defined in the region below the curve labeled qf < (?2'^^ 
in Fig. El 

In the second part of this work we have shown exam- 
ples of an alternative mechanism which is defined when 
If < 92'^^- Qualitatively, it means that the field turns 
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around on each cycle. The kinetic energy gained when 
the field is rolling down the potential in the expansion- 
ary phase, is enhanced during the collapse pushing the 
field fasrther up the potential on each cycle. The en- 
ergy density of the field comes eventually to dominate 
the curvature and background contributions establishing 
the initial conditions for inflation. 

We have seen that there are values of j for which the 
maximum value of the field is above three Planck units, 
thus able to provide a viable inflationary evolution. It 
is nonetheless observed that it is very easy to violate 
the Hubble bound Hui < 1. This problem signiflcantly 
constraints the range of parameter space for which one 
obtains a viable result within the limits of validity of the 
theory. Moreover, it involves large values of the param- 
eter j. Large values of the parameter j are considered 
to be unnatural |^. Also worth pointing out that the 
mechanism can only be implemented for potentials with 
a minimum. 

In this work we have neglected discreteness corrections. 
This approach is valid as long as d <C 1 . Admit- 

ting that the maximum value of a occurs near a*, we 
obtain that H(a = a*)a* <C 1. Using the approximation 
that the maximum of the background energy density oc- 
curs when q — I with D = I, we get the equivalent 
condition r ^ 1. From Eq. H28|l we can further extract a 
constraint on the minimum value of ginit '■ 

(y)"%.^-r^/^»l. (38) 



More specifically, ginit > 0.5 for Wc\ = 1/3, independently 
of the value of j. This bound, only involving orders of 
magnitude, is more stringent than the heuristic consis- 
tency condition used so far in the literature, Hai < 1. It 
also suggests that it is very difficult to move the field as 
far as 3mpi unless the initial kinetic energy of the field 
is at least comparable to the energy density of the back- 
ground fiuid. It would be interesting to investigate how 
the dynamics and bounds get modified once the discrete- 
ness corrections are taken into account. 

Finally we would like to keep in mind that the oscil- 
lations of the field may induce the conversion of part of 
its energy density into radiation thereby modifying the 
background's contribution in the process. The picture is 
therefore more complex than the one we dealt with here 
and deserves further study. 
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